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The electronic states at the direct band gap of monolayer WSe 2 at the K + and K~ valleys are 
related by time reversal and may be viewed as pseudo-spins. The corresponding optical interband 
transitions are governed by robust excitons. In double resonant Raman spectroscopy, we uncover 
that the 2 s exciton state energy differs from Is state energy by exactly the energy of the combination 
of several prominent phonons. Superimposed on the exciton photoluminescence (PL) we observe 
the double resonant Raman signal. This spectrally narrow peak shifts with the excitation laser 
energy as incoming photons match the 2s and outgoing photons the Is exciton transition. The 
multi-phonon resonance has important consequences: Following linearly polarized excitation of the 
2s exciton a superposition of valley states is generated which can relax fast via phonon emission 
and with minimal loss of coherence from the 2s to Is state. This explains the high degree of valley 
coherence measured for the Is exciton PL. 


Introduction .— Monolayer (ML) transition metal 
dichalcogenides (TMDCs) such as M 0 S 2 and WSe 2 have 
a direct optical bandgap at the AT-points of the Brillouin 
zone m- Light absorption and emission are governed 
by robust excitons, Coulomb bound electron-hole pairs, 
as their binding energy and oscillator strength is greatly 
enhanced in these ideal two-dimensional (2D) systems. 
The electronic states in the K + and K~ valleys are re¬ 
lated by time reversal and may be viewed as pseudo¬ 
spins. The electrons can be optically initialized either 
in the K + or K~ valley using the chiral optical selec¬ 
tion rules, i.e. with circularly right (ct + ) or left (ct - ) 
polarized light, respectively HE]. Hence, the valley de¬ 
gree of freedom can be manipulated by light Elea , an 
opportunity to investigate Berry phase effects in solid- 
state physics El- Highly circularly polarized photolumi¬ 
nescence (PL) reported for WSe 2 [H E], M 0 S 2 [ITHTi] 
and WS 2 PE] confirms efficient and robust optical valley 
polarization initialization. For coherent manipulation of 
arbitrary valley states, efficient interaction with linearly 
and eventually elliptically polarized laser light needs to 
be demonstrated. Excitation with linearly polarized light 
yields the coherent superposition of exciton pseudo-spin 
states, a process called optical alignment of excitons PD¬ 
As each exciton state in ML TMDCs is built from elec¬ 
trons and holes from the I\ + and K~ valleys, optically 
aligned exciton states are also referred to as coherent val¬ 
ley states [Si. However, the process of optical generation 
of valley coherence and its detection in linearly polarized 
PL emission has been more elusive mum]. 

In most semiconductor systems (pseudo-)spin coher¬ 
ence is lost during energy relaxation involving scattering 
events. As a consequence strictly resonant laser exci¬ 
tation is favoured for coherent control schemes pZHl9]. 
In this work we demonstrate how and why valley coher¬ 
ence can be generated efficiently capitalizing on a nrulti- 
phonon resonance that we uncover in double resonant 
Raman scattering p20M24] . We combine two remarkable 
properties of ML TMDCs in our experiments, namely (i) 
the dominant role of excitons 0IM31] combined with 


(ii) the strong interaction of excitons with phonons HD- 

m 

We show that valley coherence generation in ML 
WSe 2 is most efficient when the excitation laser is ex¬ 
actly 140 meV above the Is exciton. This has two rea¬ 
sons: First, light absorption is enhanced as states of the 
2s (or 2 p) exciton are available. Second, fast and effi¬ 
cient energy relaxation with minimal coherence loss to 
the Is exciton occurs via multi-phonon emission. We ob¬ 
serve that the exciton PL emission from coherent valley 
states is superimposed on a strong double resonant Ra¬ 
man scattering signal (not linked to valley physics) as 
the incoming and outgoing photon energies are in reso¬ 
nance with the 2s and Is exciton states, respectively. We 
are able to spectrally separate the exciton emission and 
the spectrally narrower Raman signal by systematically 
varying the laser energy. We show that the Raman signal 
provides a substantial contribution to the global emission 
and its linear polarization. Therefore care needs to be 
taken when analyzing strongly linearly polarized emis¬ 
sion from the Is exciton. Double resonances are much 
stronger than single resonances [40], making double res¬ 
onant Raman scattering a sensitive tool for investigat¬ 
ing the electronic transitions in, for example, graphene, 
graphite and carbon nanotubes [52HM]- Thus, the ob¬ 
servation of double resonant Raman scattering mediated 
by excitons in TMDCs opens up possibilities for in-depth 
studies of fine details of the energy spectrum, currently 
debated in the literature [3 ESH33] • 

Samples and Experimental Set-up The WSe 2 ML 
flakes are obtained by micro-mechanical cleavage of a 
bulk crystal on Si0 2 /Si substrates. Experiments at 
T = 4 K are carried out in a confocal microscope m- 
The detection spot diameter is ss 1 //in, i.e. considerably 
smaller than the ML size of ~ 10 /mix 10 /mi. For time 
integrated experiments, the PL emission is dispersed in a 
spectrometer and detected with a Si-CCD camera. Our 
target is to distinguish between PL emission and spec¬ 
trally narrow Raman signals. The WSe 2 ML is there¬ 
fore excited by picosecond pulses generated by a tunable 
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FIG. 1: 2-photon excitation (a) Broad 2s exciton PL 
peak following excitation of 3p state. Trion (T) and neutral 
exciton X° (Is) are marked, (b) emission at the X° energy (Is 
exciton) for different laser energies 2 El close to the 2 p exciton 
resonance. Laser polarization linear X, detection linear X (co¬ 
polarized). (c) same as (b) but for detection linear Y (cross- 
polarized). (d) same data as (b) but presented as a contour 
plot to show changes in PL intensity and shape as a function 
of 2 El, dashed line marks 2 p resonance, (e) same data as 
(c), same colour scale as (d). 


frequency-doubled optical parametric oscillator (OPO) 
synchronously pumped by a mode-locked Ti:Sa laser. 
The typical pulse temporal and spectral widths are 1.6 
ps and 3 meV, respectively, the latter being smaller than 
the exciton PL emission linewidth ~ 15 meV; the repeti¬ 
tion rate is 80 MHz. 

Results and Discussion .— In our measurements with 
a tunable laser source, we change the excitation energy 
and monitor the PL emission of the WSe 2 ML. We con¬ 
centrate here on the emission of the neutral exciton X° 
ground state (Is) at 1.75 eV, with a FWHM of about 
~ 15 meV, see Fig [!}). At 1.72 eV we observe the charged 
exciton emission (trion, T). The emission observed at 
lower energies is probably related to localized states, it 


disappears when the temperature is raised SB- In previ¬ 
ous work [9], we have identified the position of the 2s and 
the 2 p exciton states in 1-photon and 2-photon PL exci¬ 
tation (PLE) experiments, respectively. Interestingly, in 
ML WSe 2 , both states are very close in energy (around 
1.89 eV). Here we show in addition hot luminescence of 
the 2s state in Fig. [L}i as a result of resonant 2-photon 
excitation of the 3 p state at about 2.03 eV. Crucially 
for the experiments presented here, resonant excitation 
of the excited (2s/2p) exciton state results in a strong 
enhancement of the Is exciton PL. Importantly, using a 
linearly polarized laser to excite the 2s/2p state, the re¬ 
sulting Is emission is strongly linearly polarized ~ 40% 
(compare intensities in Figjlji vs. Fig[lj3 and data for 
1-photon PLE in the supplement [4(|j), a dear finger¬ 
print of efficient valley coherence generation i.e. exciton 
alignment [9]. It is very surprising to observe the recom¬ 
bination of a coherent superposition of valley states using 
optical excitation as high as 140 meV above the Is state. 

Our 2-photon PLE experiments shown in FigsjlJj-e 
clearly demonstrate that the variation of the exciton laser 
energy has several effects on photons emitted at the Is 
energy: (i) excitation close to the 2 p exciton results in 
drastic changes of the X° PL emission intensity (ii) as the 
laser is scanned through the 2 p resonance, the line shape 
and central energy of the emission are modified. Equiv¬ 
alent results for 1-photon PLE are shown in the supple¬ 
ment m • A detailed analysis of the shape and energy 
position of the X° emission for different laser energies is 
shown in Fig. for 1-photon PLE. When exciting close 
to the 2s resonance, the emission consists of two distinct 
features: In addition to the strong and highly polarized 
X° emission (energy position confirmed for more than 50 
values of non-resonant laser excitation energies) we ob¬ 
serve a spectrally much narrower feature, that changes 
energy as the laser energy is changed. Careful fitting of 
the emission for different laser energies confirms that the 
additional, narrow feature is observed exactly 140 meV 
below the laser excitation energy, as shown clearly in 
Fig. In Fig. we show that the linewidth of the 
X° PL emission is about 15 meV and the narrower con¬ 
tribution has a FWHM of about 5 meV [22]. We have 
confirmed this observation in several samples. In Fig. 
we show additional evidence for a different sample: Us¬ 
ing 2-photon excitation, we also find a strong resonance 
when twice the laser energy 2 El is 140 meV above the 
Is exciton. Again we analyse the X° emission around 
this resonance: Two features, one moving with the laser 
energy the other one not, can be discerned. Fitting the 
results in Fig. clearly reveals a narrow feature that is 
always 140 meV below 2 El, in addition to the X° PL 
with fixed emission energy, see Fig. Hi. 

Our experimental results raise the following questions: 
What is the physical origin of the spectrally narrow peak 
whose energy is fixed with respect to that of the laser? 
Why is the X° coherence and emission intensity at a max¬ 
imum when this narrow peak centre energy is exactly at 
the X° Is resonance? A maximum in the PLE signal re- 
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FIG. 2: 1ML WSe 2 , T=4 K. (a) emission at exciton X° Is energy (red circles) for different excitation laser energies El for 
sample 1, spectra are offset for clarity. Two contributions can be separated via fitting: The X H PL component (green) and 
the resonant Raman contribution (pink, shaded grey), the latter changes with El from top to bottom panel, (b) Extracted 
transition energies for X° PL (green squares) are constant, resonant Raman spectral position (open circles) changes linearly as 
a function of laser excitation energy, for comparison we plot El — 140 meV (dashed line), (c) and (d) same as (a) and (b) but 
using 2-photon excitation for sample 2. (e) Line width of the X° PL (green squares) and resonant Raman signal (pink circles) 
as a function of laser excitation energy as extracted from panel (a), (f) Simply resonant Raman scattering as outgoing photons 
are in resonance with the broad 2s exciton state shows 3 phonon replica more than 90 meV below the 2-photon excitation, (g) 
schematics of exciton generation and relaxation via multiple phonon processes, leading to enhanced X° emission and a double 
resonance (incoming and outgoing photons) for Raman scattering. 


quires (a) strong absorption by an excited state, here the 
2s or 2 p exciton, and (b) efficient relaxation from the ex¬ 
cited state to the exciton ground state (Is). Relaxation 
between electronic states in semiconductors is most effi¬ 
cient by emission of phonons. This allows in addition to 
partially preserve coherence 55j- Raman studies in ML 
TMDC have revealed a multitude of well characterized 
phonon modes, e.g. see [551155115511551 . A key point in 
our experiment is the enhancement of light scattering by 
phonons as the photon energy is in resonance with an 
electronic transition. In Fig. E we present an example 


of resonant Raman scattering in our sample using a sin¬ 
gle (not double) resonance, similar to recent experiments 
in ML M 0 S 2 [35]: Excitons are generated via 2-photon 
absorption, and up to three optical phonons (with en¬ 
ergy 31 meV [55] ) are emitted, so that excitonic final 
states fall into the broad line of the 2s/2p state. This sim¬ 
ple resonant Raman experiment shows three clear phonon 
replica more than 90 meV below the excitation energy. 
Multiple Raman features are often observed in 2D mate¬ 
rials such as GaSe [50] and TMDCs [36] as the interaction 
of phonons with exciton transitions is efficient. Due to 
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TABLE I: Selected phonon energies in WSe 2 , after Ref. [36] 


Mode 

Raman shift [cm *] 

Energy [meV] 

E'( r) 

248 

30.73 

A'lCn 

250 

30.98 

3 LA(M) 

394 

48.82 


2D confinement, exciton oscillator strength and binding 
energy are enhanced and emission linewidths in the meV 
range allow to see clear resonance effects. For 2s/2p to 
Is phonon-assisted transitions investigated here in ML 
WSe 2 an ideal situation for double resonant Raman scat¬ 
tering occurs: The incoming photons are in resonance 
with the 2s/2p exciton state and the outgoing photons 
with the Is exciton state, see Fig. [2]g. Double resonant 
Raman scattering is extremely efficient as the scattering 
cross sections are strongly enhanced nmnnsa, as doc¬ 
umented already in the pioneering experiments in GaAs 
based structures EDI EH- a model description of double 
resonant Raman processes for 1 and 2-photon excitation 
experiments is given in the supplement [HD] . 

Our next target is to analyse which phonon modes can 
be combined to bridge the 140 meV energy difference 
between the 2s/2 p and Is exciton states, respecting the 
wavevector conservation and symmetry (see supplement 
for details HO])- The WSe 2 monolayer has D 3 h point 
symmetry and P6m2 (187 or D^ h ) space group. The 
analysis of Refs. [44H46] [54l allows to associate 9 phonon 
modes at the T point (center of the Brillouin zone) 
with the following representations of D 3 h point group: 
A[ + 2A" + 2E' + E". To be Raman-active, the phonon 
modes should transform as quadratic combinations of co¬ 
ordinates, i.e. according to the irreducible representa¬ 
tions A' x , E', and E" m • Theoretical calculations of 
phonon energies in WSe 2 monolayers are presented in 
Refs. [351 113 and peaks in Raman spectra are identi¬ 
fied in Refs. J35J 135] . with strong signals with well de¬ 
fined peaks attributed to the A\ (r), E'(T) and 3LA(M) 
modes. The combinations of the prominent phonons 
A' 1 {r)®2E'(T)®3LA(M), or 2A' 1 (r)0.E , (r)0 3LA(M), 
or 2>A! X (T) 0 3LA(M), or 3 E'(T) 0 3LA(M) [53] have an 
energy of hu> sa 140 meV, see Table [I] This is very close to 
the energy separation between 2s/2p and Is excitons in 
ML WSe 2 and can account for double resonant Raman 
scattering [HDl. These combinations contain the modes 
transforming both according to A[ and E' and, hence, are 
Raman-active in z(xx)z and in z(xy)z configurations, i.e. 
both in co- and cross-linear polarizations. Analysis shows 
that these multiphonon modes are consistent with both 
resonant Raman transitions from 2s to Is and from 2 p 
to Is exciton states observed, respectively, in one-photon 
and two-photon excitation experiments [40] . 

Our detailed analysis of double resonant Raman scat¬ 
tering allows to distinguish between the different contri¬ 
butions to the ML WSe 2 emission following excitation 
with linearly polarized laser light. As we can spectrally 
separate the X° PL from the Raman signal, we can ana¬ 



FIG. 3: (a) Following linearly polarized laser excitation 

140 meV above the X° PL emission energy, the PL is strongly 
polarized, (b) The narrower double Raman resonance is also 
strongly linearly polarized. The sum of the graphs shown in 
(a) and (b) is a least square fit for the global emission with 
Gaussian line-shapes. We plot them here in different panels 
for clarity. 


lyze the polarization of each contribution individually, as 
shown in Fig.[3j For a laser energy 140 meV above the X° 
PL energy we observe a maximum of the global emission 
intensity, but also for each of the two components taken 
individually. For the X° emission, we observe a linear po¬ 
larisation degree of 42%, corresponding to strong valley 
coherence i.e. exciton alignment. From this time inte¬ 
grated measurement we can infer that the coherence is 
maintained during the PL emission time that we have de¬ 
termined to be in the ps range m- For the Raman signal, 
we observe a linear polarization degree of 30%. The de¬ 
polarization of Raman line is qualitatively in agreement 
with the selection rules for the relevant modes, see above. 
In our experiment we successfully separate PL and Ra¬ 
man contribution to the signal, an important distinction 
for valley coherence control schemes. For other semicon¬ 
ductors this separation can be achieved in time-resolved 
experiments, as the scattered light is only present during 
the laser pulse whereas the PL signal persists for longer. 
This is not possible here as the PL signal decays within 
ps, i.e. during the excitation pulse of our ps laser, here 
experiments with fs time resolution would be helpful. Al¬ 
though the same physical processes analysed in this work 
will be applicable to experiments with femto-second laser 
pulses, in the spectral domain the Raman signal will ap¬ 
pear much broader than the PL linewidth and will not 
be separable. 

Conclusions .— We analyze in detail the neutral exci¬ 
ton emission (Is) and the scattered light detected when 
exciting the 2s (or 2 p) exciton state with a linearly polar¬ 
ized laser. Tuning the laser in resonance with the 2s (or 
2 p) state has two effects: First, we observe double reso¬ 
nant Raman scattering. The incoming photons are reso¬ 
nant with the excited exciton state, the outgoing photons 
with the Is exciton state. Second, we observe strong ab¬ 
sorption, followed by efficient energy relaxation via emis¬ 
sion of several phonons, as for example the combination 
Ai(r) 0 2E'{T) 0 3LA(M). As a result the Is exciton 
emission is strongly enhanced as compared to excitation 
with other laser energies. The PL emission is linearly 
polarized above 40%, corresponding to very efficient val- 
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ley coherence generation following the resonant exciton 
generation. Importantly, we are able to clearly separate 
the PL emission from the spectrally narrow Raman sig¬ 
nal as we vary the laser excitation energy, which shifts 
the Raman feature.The multi-phonon resonance between 
the excited exciton states might be one of the reasons 
why efficient valley coherence can be generated in ML 
WSe 2 , contrary to M 0 S 2 or MoSe 2 , where this resonance 
might not occur. In particular we do not observe the 
double resonant Raman feature in PLE experiments on 
high quality MoSe 2 MLs around the 2 p resonance f55| . 
Further theoretical work is needed to reveal the nature 
(perturbative or strongly coupled) of the exciton-phonon 
interaction 153 ■ Also the influence of magnetic fields 


applied perpendicular to the sample plane will allow to 
study separately the evolution of the polarization of the 
Raman signal and the valley coherence of the neutral ex¬ 
citon [551 152] . The exact energy spacing and symmetry 
of the exciton states that govern the optical properties of 
ML TMDCs are still under debate with scattered values 
in the literature [9J H5H3U • The demonstration of double 
resonant Raman scattering mediated by excitons intro¬ 
duces a sensitive tool to probing this material. 
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SUPPLEMENTARY INFORMATION 

I. ADDITIONAL EXPERIMENTAL DATA: ONE 
PHOTON PHOTOLUMINESCENCE SPECTRA 


Figure [4] shows the 1-photon photoluminescence (PL) 
spectra measured at the conditions similar to those pre¬ 
sented in the main text for two-photon PL. For these 
measurements the excitation laser energy was scanned in 
the vicinity of 2s exciton energy, 1.89 eV. The change of 
shape of the PL emission as well as the variation of the 
position of the central energy with the energy of excita¬ 
tion laser is clearly seen. 


Laser Lin. X /Det. Lin. X Laser Lin. X /Det. Lin. Y 
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FIG. 4: 1-photon excitation (a) emission at the X° energy 
(Is exciton) for different laser energies El close to the 2s ex¬ 
citon resonance. Laser polarization linear X, detection linear 
X (co-polarized). (b) same as (a) but for detection linear Y 
(cross-polarized), (c) same data as (a) but presented as a 
contour plot to show change in PL intensity and shape as a 
function of the excitation laser energy El- (d) same data as 
(b), same color scale as (c). Polarization configurations are 
marked on top of the panels. 


II. SYMMETRY OF PHONON MODES IN WSE 2 
AND SELECTION RULES FOR RAMAN 
PROCESS 


WSe 2 monolayer has D point symmetry and P6m2 
(187 or D\ h ) space group pQ. There are three atoms 
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in the unit cell resulting in 9 phonon modes. The de¬ 
tailed symmetry analysis of the modes is presented in 
Refs. [IHH]. The symmetry of the phonon modes at the 
center of Brilloin zone can be established using correla¬ 
tion method. We choose a center of point transforma¬ 
tions in the position of W atom. The symmetry of W 
site is, hence, D 3hl while the symmetry of the Se sites is 
C 3v , see Fig.[5ja). In D 3 h point group the displacements 
transform according to E' representation (x , y) and A 3 
representation (z), see Table |TT} In C 3v point group the 
displacements transform according to E (x, y) and A 3 
(z). Making use of the correlation table for D 3 h point 
group [B] we obtain the following representations of the 
phonon modes mmmm- 

A 3 +E’ + E’ + E" + A[ + A '2 = 2A" + 2E' + E" + A[. (1) 

The acoustic modes at the T point transform according 
to Aij (non-degenerate) and E' (two-fold degenerate at 
the T point), the remaining modes are optical. 



FIG. 5: (a) Lattice of WSe 2 projected on the Tungsten plane, 
(b) Brilloin zone of WSe 2 . 

We recall that the Raman scattering can be described 
as phonon-induced fluctuations of polarizability tensor 
X which links the linear dielectric polarization of the 
medium P and electric field E as 

Pi = XijPji (2) 

where i and j denote Cartesian components. For a 
single-phonon process the fluctuating part of Sxij is 
proportional to exp (=piflQ,t), where a enumerates 
phonon modes (including their polarization and wavevec- 
tor), are the generalized phonon coordinates, L l a are 
their frequencies, and =F signs correspond to the absorp¬ 
tion/emission of a phonon, respectively. Analogously, for 
multiphonon processes involving a, f3,... ,5 modes, we 
have 

S\ij oc Q (a) Q (/3) ... exp [-i(±O a ±Clp±...± f2j)i], 

( 3 ) 

giving rise to the combination frequencies in Raman spec¬ 
tra. To be Raman-active, the product Q^QW ... Ql* 5 ) 


should transform like products of coordinates r ; ;'r 7 , where 
r = (; x,y,z ). Hence, the Raman-active modes in WSe 2 
transform as quadratic combinations of coordinates, i.e. 
according A[, E 1 , and E". Hence, the only Raman- 
inactive mode is A% but the combined modes, e.g., 2 A' 2 ' 
are also active (A' 2 ' x A 2 = A \). Note, that the modes 
with E" symmetry are Ranran-active at the oblique in¬ 
cidence only. 

TABLE II: Representations of D 3 h point symmetry group 


Representation 

Basic functions 

A'i 

x 2 + y A ; z 2 

A^ 

Jz 

E' 

{x,y)', ( x 2 — y 2 , —2xy) 

A'i 

zJ z 

AS 

z 

E" 

( xz,yz) 


We are also interested in the symmetry of the phonons 
at the M point (K || y) of the Brillouin zone, Fig. [5](b) . 
The wavevector group is C 2 V {&h becomes a v and one 
of a v planes, namely, containing T-M or y direction, re¬ 
mains). In C' 2 v point group displacements transform ac¬ 
cording to A\ (invariant, y), B\ (z) and B 2 ( x ) represen¬ 
tations. Particularly, the longitudinal acoustic mode at 
the M-point, LA{ M), transforms according to A\ in C 2 V 
(note, that E' = A 3 + B 2 ). 

The multiphonon modes A^r) (g) 2£ I, (F) g) 3LA(M), 
or 2Ai(r) 0 E'(T) 0 3LA(M), or 3Ai(F) 0 3LA(M), 
or 3£’ , (r) 0 3LA(M), involve combinations of optical 
phonons A^r), E'(T) at the Brillouin zone center and 
a combination of three longitudinal acoustic modes at 
M-point with zero net wavevector. They contain several 
A[ and E' representations. Hence, these combinations 
are Raman-active in z(xx)z configuration and in z(xy)z 
configurations. Note, that due to the fact that these mul¬ 
tiphonon states contain A( and E' representations they 
can result in the transfer of excitation from 2 p exciton 
state (representation E' x E' = A\ + A' 2 + E', where the 
first representation in the product denotes the symmetry 
of the exciton Bloch function and the second one denotes 
the symmetry of the envelope) to Is exciton state (rep¬ 
resentation E' x A' = E'). 

III. MODEL OF DOUBLE RESONANT RAMAN 
SCATTERING 

To illustrate the enhancement of the Raman scattering 
we model 2s and Is excitonic states in WSe 2 as in-plane 
dipoles with the polarization Pis and P\ s obeying har¬ 
monic oscillator equations [8] 

P 2s + ojj s P 2s = gE 0 cos ut, (4a) 

Pis + u\ s Pis = S(t)P 2s , (4b) 

where a> 2 s and Wi s are the energies of the correspond¬ 
ing excitonic states, Eq cos cot is the electric field of the 


















incident radiation, whose frequency ui is assumed to be 
close to 2 s exciton resonance, g is the coupling strength 
and the matrix S(t) describes multi-phonon-induced in¬ 
teraction of the oscillators. Its matrix elements in Carte¬ 
sian basis are S'y = sfj' 1 cos fit, where f 2 is the frequency 


of the multi-phonon mode [ 12 ] and the prefactors 
contain contributions from virtual transitions involving 
all relevant phonon modes. Their calculation is highly 
complex since it involves high order perturbation theory 
and is beyond the scope of the present work (13]. Equa¬ 
tion (4b) clearly shows that only (multi-)phonon modes, 
which transform as products of coordinates, x 2 , y 2 , xy 
or their combinations (that are A! x and E’) can couple 2s 
and Is excitons. 

The solution of set of Eqs. © yields the polarization of 
the Is exciton. It contains the contributions oscillating 
as cos [(w ± 12 ) 2 ] with the amplitudes 


„ SffgEoj cos [(w ± 12 ) 2 ] ^ 

t^«- w2 )[<-( w±n ) 2 ]' 

Equation ([5]) shows that the polarization of Is exciton os¬ 
cillates with the frequencies ui ± 12 which corresponds to 
the anti-Stokes and Stokes peaks [14]. The denominators 
in Eq. 0 clearly demonstrate the resonant enhancement 
of the Raman scattering m The double resonance oc¬ 
curs if incident radiation frequency is close to ui 2 s and, 
additionally, if ui i s is close to w —12 (or w + O for the anti- 
Stokes scattering) and results in drastic enhancement of 
Pis and, correspondingly, of the scattered light intensity. 

Similar analysis demonstrates that the double reso¬ 
nance Raman scattering is possible under the conditions 
of two-photon excitation, where the 2 p exciton resonance 
is addressed. Note, that in WSe 2 the 2s and 2 p excitons 
are very close in energy [9| . The 2 p exciton state can be 
conveniently modelled in the coupled-oscillator picture as 
a quadrupolar-like mode U 2 Pl ij, where one of the Carte¬ 
sian subscripts (2) denotes the Bloch function of exciton, 
representation E’ , and the other one (j) denotes the en¬ 
velope function (p^-shell or p y -shell like, representation 
E ') for 2 p state. The set of Eqs. |4]) is replaced by 

P 2 p,ij T = 'y fJi j,ki Co./ cos . ( 6 a) 

kl 

As,j + W l s Pls,i = ^ S'i,jl(t)U2p,jl, ( 6 b) 

jl 


where the tensor g C describes the selection rules at two 
photon absorption and components S[ ^ oc cos 122 de¬ 
scribe contributions for multiphonon processes. For in¬ 
plane oscillators where 2, j, l = x or y the phonons of E' 
and A! x symmetry can provide such a coupling, i.e. the 
same modes which provide 2s to Is transfer. The solu¬ 
tion of Eqs. 0 can be presented in a form, similar to 


Eq. j5j: 


p ^iJl9jl,km^O,kEo,m COS ](2w ± 12)2] 

K-4^)[<-(2. ± I2) 2 ] ; U 
Here the double resonance occurs when the twice the 
photon energy 2 w approaches W 2 P , the frequency of the 
2 p state, and the frequency of the secondary light, 2 w±f 2 , 
approaches the frequency of the ground exciton state, Is. 

The experimental data presented in the main text and 
in Fig. [4] clearly show that the double resonance Raman 
signal is superimposed on the PL. To analyse the origin of 
these two contributions to the emission and demonstrate 
that such a superposition is indeed consistent with the 
model presented above, let us consider the correlation 


function of the source terms for P ls in Eqs. (4b), ( 6 b). To 


be specific, we consider the case of 1 -photon excitation, 
Eqs. Q, and introduce the correlation functions as 


= Jl (cos [u 2 S (t + T)]S ik (t + r)P k (t + T)Sji(t)Pi(t)), 

kl 

( 8 ) 

where the averaging is performed over the time 2 at a 
fixed delay r; the factor cos ( uj 2 S t ) is included into Eq. ^ 
to exclude the oscillations of S^-(t) at the optical fre¬ 
quency. The correlation functions 3^ quantify the ex¬ 
tent to which the driving force in the right hand side 
of Eq. (4b I is harmonic. Particularly, 3jj can be pre¬ 


sented as a sum of contribution of all possible rnulti- 
phonon modes. The terms of interest for us oscillate at 
a frequency 12 ~ u> 2 S — wi s , which provide the double 
resonance effect: 


j(r) = 


K —0 


cos (fir) exp (— 7 qt). 


(9) 


Here the subscript a runs through all relevant multi¬ 
phonon combinations, see Sec. [Hj and we introduced the 
mode damping rates 7 a , small differences of the multi¬ 
phonon combination frequencies are disregarded, possi¬ 
ble effect of phonon modes dispersion is included into 
the damping rates 7 a . 

For the modes with H/ 7 a 1 the correlation function 
oscillates with time and weakly decays. These modes pro¬ 
vide double resonant Raman scattering as the induced 
polarization P ls oscillates with time as cos [(w 2 p — 12 ) 2 ] 
(or as cos ](w 2 P + L2)t] for anti-Stokes scattering). Other 
modes, where 1 $ 1 , do not act as a harmonic driving 
force for the ls-exciton polarization, since their correla¬ 
tion function decays at a time scale, 7 ” 1 , being short as 
compared with the oscillation period. These modes act 
as random (Langevin) forces and result in the incoher¬ 
ent population of Is state, giving rise to oscillations of 
Pis with the Is exciton frequency Wi s and yielding the 
photoluminescence superimposed on the Raman signal. 
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